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Abstract 

We discuss a new estimate of e'/e in the kaon system. The present approach is based on the 
evaluation of the hadronic matrix elements of the AS = 1 effective quark lagrangian by means 
of the chiral quark model, with the inclusion of meson one-loop renormalization and NLO 
Wilson coefficients. The estimate here reviewed is fully consistent with the AI = 1/2 selection 
rule in K —>■ tttt decays, which is well reproduced within the same framework. 



''Talk given at the Workshop on K Physics, Orsay, France, May 30 - June 4, 1996. 



1 Introduction 



The real part of e'/e measures direct CP violation in the decays of a neutral kaon in two 
pions. It is a fundamental quantity which has attracted a great deal of theoretical as well as 
experimental work. Its determination answers the question of whether CP violation is present 
only in the mass matrix of the neutral kaons (the superweak scenario) or is also at work directly 
in the decays. 

On the experimental front, the present results of CERN (NA31) and Fermilab (E731) are 
tantalizing insofar as the superweak scenario cannot be excluded and the disagreement between 
the two outcomes leaves still a large uncertainty. 

On the theoretical side, much has been accomplished, although the intrinsic difficulty of a 
problem that encompasses scales as different as m t and m n weights against a decisive progress 
in the field. 

The short distance (perturbative) QCD analysis has greatly profited in recent years from 
the work done by the Munich IJ and Rome groups who computed the anomalous dimension 
matrix of the ten relevant operators to the next-to-leading order (NLO), thus reducing the 
uncertainty related to a s at the 10% level. 

The largest uncertainty arises in the long-distance part of the effective lagrangian, the 
computation of which involves the evaluation of the hadronic matrix elements of the quark 
operators. 

There exist two complete estimates of such hadronic matrix elements performed by the 
aforementioned groups, and recently updated in ref. || for the lattice (at least for some of the 
operators). A recent review of 1/N C results is given in ref. 

Both groups seem to agree on the difficulty of accommodating within the standard model 
a value substantially larger than 10 -3 and obtain 10~ 4 as the preferred scale for e'/e. This 
unexpectedly small value is the result of the cancellation between gluonic and electroweak 
penguin operators ||. This cancellation, which depletes by an order of magnitude the natural 
magnitude oie'/e in the SM, marres any theoretical attempt of predicting direct CP violation 
with a precision better than a factor two (due solely to the intrinsic short- distance uncertainty). 

In view of that, it seems to us that independent estimates of e'/e, even in phenomenological 
models, are desirable. In our opinion it is crucial that a reliable evaluation of the hadronic 
matrix elements first provides a consistent picture of kaon physics, starting from the CP- 
conserving amplitudes and, in particular, by reproducing the AI =1/2 selection rule, which 
governs most of these amplitudes as well as the quantity e'/e itself. We also feel that the 
same evaluation should pay particular attention to the problem of achieving a satisfactory 75- 
scheme and scale independence in the matching between the matrix elements and the Wilson 
coefficients, the absence of which undermines any estimate. 

In ref. ||, hereafter referred as I, we have completed the study of the hadronic matrix 
elements of all the ten operator of the AS* = 1 effective quark langrangian by means of the 
chiral quark model (xQM) |7j and shown in ||, hereafter referred as II, that the inclusion of 
non-perturbative 0(a s N c ) corrections and the SU(3) breaking effects in the one- loop meson 



renormalization provided an improved scale independence and, more importantly, a good fit of 
the AI = 1/2 selection rule (see also the contribution of M. Fabbrichesi at this workshop). 

In a recent paper || we presented a new computation of e'/e based on the xQM approach. 
Our estimate takes advantage, as the existing ones, of 

• NLO results for the Wilson coefficients; 

• up-to-date analysis of the constraints on the mixing coefficient Im X t . 
Among the new elements introduced, the most relevant are 

• A consistent evaluation of all hadronic matrix elements of the ten effective quark operators 
in the xQM (including next-to- leading O(N) and 0(a s N) contributions). 

• Inclusion in the AS" = 1 chiral lagrangian of the complete bosonization 0(p 2 ) of the 
electroweak operators Qj and Q$. Some relevant 0(p 2 ) terms have been neglected in all 
previous estimates; 

• Inclusion of the meson-loop renormalization and scale dependence of the matrix elements; 

• Consistency with the AI = 1/2 selection rule in kaon decays; 

• Matching-scale and 75-scheme dependence of the results below the 20% level. 

Our prediction of e'/e depends most sensitively on the value of the quark condensate, the 
input parameter that controls penguin- diagram physics. For this reason, we discuss a inclusive 
estimate based on a conservative range of (qq), as well as the variations of all the other inputs: 
m t , lm\ t (which depends, beside m t and m c , on Bk and other mixing angles) and Aqcd- Such 
a procedure provides us with the range of values for e' / e that we consider to be the unbiased 
theoretical prediction of the standard model. Unfortunately, this range turns out to be rather 
large, spanning from —5 x 1CT 3 to 1.4 x 10~ 3 . On the other hand, it is as small as we can get 
without making some further assumptions on the input parameters, assumptions that all the 
other available estimates must make as well. 

Such uncertainty notwithstanding, we agree in the end with refs. |j] that it is difficult to 
accommodate within the standard model a value of e'/e larger than 10~ 3 . As a matter of fact, 
our analysis points to definitely smaller values, when not negative. This can be understood 
not so much as a peculiar feature of the xQM prediction as the neglect in other estimates of a 
class of contributions in the vacuum saturation approximation (VSA) of the matrix elements of 
the leading electroweak operators, which enhance the destructive interference between gluonic 



and electroweak penguins. This problem is discussed in detail in I and in ref. p0| . These new 
contributions are responsible for the onset of the superweak regime for values of m t smaller 
than 200 GeV. In our computation, it is the meson loop renormalization that in the end brings 
back e'/e to positive values. 



2 Hadronic Matrix Elements and B { Factors 



The quark effective lagrangian at a scale p, < m c can be written as 

£ AS =i = ~^ v ud Vus E + m(j4\Qt(j*) = -^E c iWftW • (2-i) 

The Qj are four-quark operators obtained by integrating out in the standard model the 
vector bosons and the heavy quarks t, b and c. A convenient and by now standard basis 
includes ten quark operators. The functions Zi(p) and yi(p) are the Wilson coefficients and 
Vij the Kobayashi-Maskawa (KM) matrix elements; r = —VtdV* s /V u dV* s . Following the usual 
parametrization of the KM matrix, in order to determine e'/e, we only need the yi(fi), which 
control the CP- violating part of the amplitudes. 

In paper I we have computed all hadronic matrix elements of the ten effective quark oper- 
ators in eq. ( |2.1|) in the framework of the xQM. The matrix elements of the first six operators 
were first computed in the xQM in ref. |lT| . The matrix elements are obtained by the inte- 
gration of the constituent quarks by means of dimensional regularization. The loop integration 
leads to results that depend on the scheme employed to deal with 75 but are scale independent. 
The renormalization-scale dependence is introduced in our approach by the meson-loop renor- 
malization of the amplitudes, as explained in I. The meson-loop corrections together with the 
gluon-condensate contributions are the most relevant ingredients in reproducing the AI = 1/2 
selection rule in K — > nn decays, (as discussed in II). 

The xQM results are expressed in a double power expansion on M 2 /A 2 C and p 2 /A 2 c , where 
M is a dimensionful parameter of the model which can be interpreted as the constituent quark 
mass in mesons, p is a typical external momentum, and A x ~ 1 GeV is the chiral symmetry 
breaking scale. 

The value of M is constrained [jn| by experimental data on the decay of tt° and rj to be 
M = 223 ± 9 MeV (243 ± 9 MeV if higher order corrections are included). Comparable values 
are found by using vector-meson-dominance, or fitting all input parameteres in the extended 



Nambu-Jona-Lasinio model 13 



The integration of the fermion degrees of freedom leads naturally to an effective bosonic 
representation of the AS = 1 quark lagrangian with the desired chiral properties. In I we 
have constructed the complete 0(p 2 ) chiral representation of the lagrangian in eq. (|2.1| ), where 
the local quark operators Qi are represented by a linear combination of bosonic operators B a , 
namely Qi — > J2 a G a (Qi)B a . The effective quark lagrangian is therefore replaced by 

£f =1 = -% £ CM G a {Qi) B a . (2.2) 

As discussed at lenght in I, the chiral coefficients G a determined via the xQM approach are 
75-scheme dependent. While the 75-scheme dependence arises in the xQM from the integration 
of the chiral fermions, the explicit /i-dependence is entirely due to the chiral loop renormalization 
of the matrix elements: 

(b\Cf =1 \a) = £ GQisd) G a (Qi) (b\B a \a) (fx LD ) , (2.3) 



where we have labeled by a and b the initial and final bosonic states. We remark that in our 
approach the //-dependence of the chiral loops is not cancelled by higher order counterterms, 
as it is usually required in the strong chiral lagrangian. 

The renormalization scale dependence is therefore determined order by order in the energy 
expansion of the chiral lagrangian. In this respect there is no direct counterpart to the ex- 
pansion in strong and electromagnetic couplings on which the short-distance analysis is based 
and, accordingly, we refer to the explicit //-dependence in the matrix elements as to the long- 
distance (LD) or "non-perturbative" scale dependence. A purely perturbative renormalization 
scale dependence is introduced in the matrix elements by the NLO running of the quark con- 
densate, which we include whenever a comparison between values at different scales is required. 
Otherwise, quark and gluon condensates are considered in our approach as phenomenological 
parameters. 

Physical observables must not depend on the chosen scheme for handling 75 in dimensional 
regularization and the value of the renormalization scale /a. Our aim is to test whether the 
estimate of observables is consistently improved by matching the "long-distance" 75-scheme 
and /1 dependences so obtained with those present in the short-distance analysis (in particular 
we identify hsd with \xld)- Whether and to what extent such an improvement is reproduced 
for many observables and for a consistent set of parameters, might tell us how well low-energy 
QCD is modelled in the xQM-chiral lagrangian approach that we have devised. 

In II, we have shown that minimizing the 75-scheme dependence of the physical isospin 
1 = and 2 amplitudes determines a range for the parameter M between 160 and 220 MeV. In 
II, it was also found that the \x dependence induced by the Wilson coefficients is substantially 
reduced by that of the hadronic matrix elements. 

These issues become crucial for e'/e where the 75-scheme dependence induced by the Wilson 
coefficients determines an uncertainty as large as 80% when using the 1/N C hadronic matrix 
elements (see for instance ref. Jl|]) which are scheme independent. 

In order to test the \i independence of e'/e we vary the matching scale between 0.8 and 1.0 
GeV, the highest energy up to which we trust the chiral loop corrections computed in I. We 
find that, in spite of the fact that some of the Wilson coefficients vary in this range by up to 
50%, the matching with our matrix elements reduces the /i-dependence in e'/e below 20% in 
most of the parameter space. We consider this improved stability a success of the approach. 

In order to discuss our results for the hadronic matrix elements it is convenient to introduce 
the effective factors I?- ' 2 - 1 = (Qj)o,^ M / (Qi)of A > which give the ratios between our hadronic 
matrix elements and those of the VSA. They are a useful way of comparing different evaluations. 

In table 1 we collect the Bi factors for the ten operators in the isospin and 2 channels. 
The values of the Bi depend on the scale at which the matrix elements are evaluated, the input 
parameters and M; moreover, in the xQM they depend on the 75-scheme employed. We have 
given in table 3 a representative example of their values and variations. 

The values of b[ ^ and B2 show that the corresponding hadronic matrix elements in the 

xQM are, once non-factorizable contributions and meson renormalization have been included, 

(2) 

respectively about ten and three times larger than their VSA values. At the same time, B\~ 





HV 


NDR 


H = 0.8 GeV 


jti = 1.0 GeV 


fx = 0.8 GeV 


// = 1.0 GeV 


o(0) 
u l 


in k 
1U.D 


111 
11.1 


1U.D 


111 
11.1 


B 2 


2.8 


3.0 


2.8 


3.0 


(2) 

B\> 


0.52 


0.55 


0.52 


0.55 


n 2 


0.52 


0.55 


0.52 


0.55 


I ) 


9 Q 


— O.U 


— 0. / 


"i Q 
— O.U 


-04 


1.8 


1.9 


1.0 


1.1 


B^ — Ba 

u 


1.7 -=-0.61 


1.8 -7- 0.64 


1.0-7-0.38 


1.1 -7- 0.40 


7 


3.0 2.2 


3.3 -=-2.4 


2.9 -=-2.2 


3.2 -r-2.3 


8 


3.3 2.2 


3.6 -=-2.4 


3.2^2.2 


3.5 ^-2.4 


9 


3.9 


4.0 


3.5 


3.6 


1U 


4.4 


4.7 


5.6 


5.9 


i?W 


2.7 -=-1.5 


3.0 -fr 1.5 


2.7-M.4 


2.9 -T- 1.5 


-°8 


2.1 -7- 1.4 


2.3 -7- 1.5 


2.1 -T- 1.4 


2.3 -7- 1.5 


-°9 


0.52 


0.55 


0.52 


0.55 


fl (2) 
"10 


0.52 


0.55 


0.52 


0.55 



Table 1: B{ factors in the xQM (including meson- loop renormalizations) at two different scales: 
fj, = 0.8 and 1.0 GeV and in the two 75-schemes. The results shown are given for M = 220 
MeV and (a s GG) /ir — (376 MeV) 4 , while the range given for -85-8 corresponds to varying the 
quark condensate between (—200 MeV) 3 and (—280 MeV) 3 . 

(2\ 

and B\ turn out to be at most half of what found in the VSA. These features make it possible 
for the selection rule to be reproduced in the xQM, as extensively discussed in II. 



For comparison, in the 1/N C approach of ref. [I4~|, the inclusion of meson- loop renormaliza- 
tion through a cutoff regularization, leads, at the scale of 1 GeV, to = 5.2, B^ = 2.2 and 

(2} (2) 

B\ = B 2 = 0.55, a result that is not sufficient to reproduce the AI = 1/2 rule. The similar- 

(2) (2) 

ity of the values B\ = B\ = 0.55 obtained in the xQM with the corresponding 1/N C results 
is remarkable, and yet a numerical coincidence, since the suppression originates from gluon 
condensate corrections in the xQM, whereas it is the effect of the meson loop renormalization 
(regularized via explicit cut-off) in the analysis of ref. ||14|| . 

The values of the penguin matrix elements (Q3) and (Q4) in the xQM lead to rather large 
Bi factors. In the case of Q3, the xQM result has the opposite sign of the VSA result and B3 
is negative. This is the effect of the large non-perturbative gluon correction. 

Regarding the gluon penguin operator Qq (and Q5), we find that the xQM gives a result 
consistent with the VSA (and the 1/N C approach), B e (B 5 ) being approximately equal to 
two for small values of the quark condensate and one-half at larger values. It is the quadratic 
dependence (to be contrasted to the linear dependence in the xQM) °f the VSA matrix element 
for the penguin operators that it responsible for the different weight of these operators at 
different values of the quark condensate. The lattice estimate at fi = 2 GeV for these operators 
gives B 5 = B 6 = 1.0 ±0.2 B. 



(2) 

The electroweak Bi factors are all larger in the xQM than in the VSA, except for B g 10 that 
are about 1/2 in the HV and in the NDR schemes. The lattice estimate at [l = 2 GeV in this 

(2) (2) 

case yields B\% = 1.0 ± 0.2 and B 9 = 0.62 ±0.10 ||. For new preliminary lattice results see 
the contribution of G. Kilcup at this workshop. 

(2) 

The most relevant result for e'/e is the value of B% which ranges from 1.5 to 2 times that 
of the VSA. This increase is due to two independent reasons. On the one hand, we found two 
new terms in the chiral lagrangian that have not been included so far in the VSA estimate 
of the Qj t s matrix elements. The chiral coefficients of these terms are computed in the xQM 
approach — as discussed in detail in I. 

From this point of view, what we have referred to as VSA — and used in table 1 as normal- 
ization for the £^7,8 operators — is not the complete VSA result. The inclusion of the new terms 
amounts up to a 60% increase of Bfl for small values of (qq) in the chosen range and down to 
about 10% for large values; smaller effects are found in the case of -B 7 °g. On the other hand, 
the meson-loop renormalization associated with the new chiral terms is large (see I) and adds 
up to reproduce the results shown in table 1. The increase in importance of the operator Q 8 
with respect to Qq turns into a more effective cancellation between the two operators for large 
values of the quark condensate while at smaller values the gluon penguin contribution prevails. 

3 Anatomy of e'/e in the xQM 
The quantity e'/e can be written as 

(3.1) 
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where, referring to the AS = 1 quark lagrangian of eq. ( |2.1| ), 

n o = YlVi (Qi)o > n 2 = ^i/i(Qi)j + a; ^Vi{Qi)Q ^ v +v' > ( 3 - 2 ) 

i i i 

and Im X t = Im V td V t * . 

The quantity Vt v+V i includes the effect of the isospin-breaking mixing between 7r° and the 
etas. 

A range for Im X t is determined from the experimental value of e as a function of m t and 
the other relevant parameters involved in the theoretical estimate. We will use the most recent 
NLO results for the QCD correction factors ^1,2,3 (see the contribution of U. Nierste at this 
workshop) and vary the AS = 2 hadronic parameter B^ in the rather conservative range 
Bk = 0.55 ± 0.25 , that encompasses both the xQM model prediction [HJ and other current 
determinations |16| . 



This procedure gives two possible ranges for Im \ t ~ "?7 1 1 1 V^t, | 2 , which correspond to having 
the KM phase in the I or II quadrant (p positive or negative, respectively). For example, for 
m P° lc = 180 GeV (m t (m w ) ~ 183 GeV) and A^ D = 350 MeV we find 

1.1 x 10~ 4 < ImA t < 1.9 x 10" 4 and 0.75 x 10~ 4 < Im X t < 1.9 x 10~ 4 (3.3) 
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Figure 1: Histograms of the partial contributions to e'/e of the height relevant operators for 



1.3 x 1(T 4 and A§£ D 



350 MeV. 



(qq) (0.8 GeV) = (-200 MeV) 3 , mf 5 = 180 GeV, lm\ t 
Gray (black) histograms represent the contribution of each operator without (with) meson-loop 
renormalization. The last two histograms correspond to the sum of all contributions. 



in the first and second quadrant respectively. For the range of Bk considered, varying all the 
other parameters (including m t and Aqcd) affects the above limits on Im X t by less than 20%. 
In particular, the upper bound on lm\ t is stable, becoming a sensitive function of the input 



parameters only if we consider Bk > 0.5. In other words, we agree with ref. |hJ that it is the 
theoretical uncertainty on the hadronic AS = 2 matrix element that controls the uncertainty 
on the determination of ImA*. 

It is useful to consider the individual contribution to e'/e of each of the quark operators. 
We have depicted them as histograms, where the grey (black) one stands for the contribution 
before (after) meson-loop renormalization. Henceforth all results are given for M = 220 MeV 
in the HV scheme. 

It is clear from the histograms of Fig. 1, 2 and 3 that the two dominating operators are 
Qq and Q$. Yet, since they give contributions approximately of the same size and opposite in 
sign, the final value turns out to be relatively small and of size comparable to that of most of 
the other operators. This result is at the origin the large theoretical uncertainty as well as the 
unexpected smallness of e'/e. 

The same histograms serve the purpose of showing that the meson-loop renormalizations 
are crucial not only in the overall size of each contribution but also in determining the sign of 
the final result (see Fig. 2). These corrections are here consistently included in the estimate 
for the first time. 

The role of the operator Q 4 turns out to be marginal in our approach. In comparing 
this result with that of the 1/N C framework Q (see also the final tables in ref. [|l7j where 
we reproduce the individual 1/N C contributions for the standard ten operators), it should be 
recalled that in the above analysis the Q4 operator is written in terms of Qi, Q2 and Q3 and 
that its values is therefore influenced by the B^ factors assigned to the former matrix elements. 
In particular, while B\ and B 2 are in ref. M requested to be large in order to account for 
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Figure 2: Same as Fig. 1 for (qq) (0.8 GeV) = (-240 MeV) 3 . 
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Figure 3: Same as Fig. 1 for (qq) (0.8 GeV) = (-280 MeV) 3 . 

the AI = 1/2 rule, B 3 is assigned the value of 1. Such a procedure produces a rather large 
value for the matrix element of Q4. In our approach, we see that in fact also B% is large (and 
negative!) and that Q4, once written in terms of the other operators, is small, as found in the 
direct estimate. 

We hope that this discussion has convinced the reader that the quantity e'/e is difficult 
to estimate even with a factor two uncertainty. We think that only the order of magnitude 
can be predicted in a reliable manner. The reason is very simple: the final value is the result 
of the cancellation between two, approximately equal in size, contributions. Accordingly, any 
uncertainty is most likely amplified by an order of magnitude and we are unfortunately dealing 
with rather large ones. It is also important to understand that these considerations hold for 
any theoretical approach. 

Keeping this in mind, by varying all parameters in the allowed ranges and, in particular, 
taking the quark condensate, which is the major source of uncertainty, between (—200 MeV) 3 





and (-280 MeV) 3 we find 

-27xl(T 4 <e'/e< 9xl0~ 4 , (3.4) 
where we have kept Aqq D fixed at its central value. A larger range, 

- 50 x 10" 4 < e'/e < 14 x 10" 4 , (3.5) 

is obtained by varying Aqq D as well. 

It should be stressed that the large range of negative values that we obtain is a consequence 
of two characteristic features of the matrix elements: i) the enhancement of the size of the 
electroweak matrix elements (Qsj) due to the coherent effects of the additional 0(p 2 ) contri- 
butions so far neglected (see discussion in sect. 5) and the chiral loop corrections; ii) the linear 
dependence on (qq) of the leading gluon penguin matrix elements compared to the quadratic 
dependence of the leading terms in the electroweak matrix elements, which makes the latter 
prevail for large values of the quark condensate. The effect of i) represents an enhancement of 
the leading electroweak matrix elements up to a factor two with respect to the vacuum insertion 
approximation and present 1/N C estimates (see table 3), while feature ii) is absent in the 1/N C 
approach, the quark condensate dependence being always quadratic. 

To provide a somewhat more restrictive estimate we may assume for the quark condensate 
the QCD Sum Rules improved PCAC result @, namely (qq) = -(221 ± 17 MeV) 3 at our 
matching scale \x = 0.8 GeV, and thus find 

4.5 ±ji x 10- 4 quadrant I 
3.9 it* x 10~ 4 quadrant II. 

whose average gives e'/e — (4 ± 5) X 10~ 4 . 

The range of the quark condensate on which the above estimate is based, is not favorite 
by our analysis of the AI =1/2 selection rule in the xQM. Higher values accommodate more 
naturally the rule, at least for a constituent mass M ~ 220 MeV — the value at which we also 
find 75-scheme independence of e'/ e. For values of the quark condensate in the range of Figs. 2 
and 3 the central value of e'/e shifts toward the superweak regime, and the role of meson loop 
corrections becomes crucial. For instance, by taking the simple PCAC result for the strange 
quark condensate, (qq) = —(244 ± 9 MeV) 3 at \x = 0.8 GeV, we find 

1.4 x 10~ 4 quadrant I 
1.2 ±£0 x 10- 4 quadrant II. 

Actually, for such a range of (qq), negative central values oie'/e in both quadrants are obtained 
due to the extra terms of the bosonization of the electroweak operators Qr and Qg neglected in 
the previous estimates. Only after the inclusion of the meson-loop renormalization e'/ e turns 
to the positive central values of eq. (|3.7|). 

To have an overview of the present theoretical status for e'/e in Fig. 4 we have summa- 
rized the predictions of three theoretical approaches and compared them with the present 1 a 
experimental results (labeled bands). 

The smaller error in the lattice estimate originates in the Gaussian treatment of the uncer- 
tainty in the input parameters with respect to the flat la error included in the other estimates. 
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Figure 4: Present status of theoretical predictions and experimental values for e'/e (in units 
of 10~ 4 ). The most recent 1/N C [4] and lattice [3] estimates are compared to (a) our unbiased 
estimate (3.5), (b) our more restrictive estimate (3.6). 

4 Outlook 

Our phenomenological analysis, based on the simplest implementation of the xQM an d 
chiral lagrangian methods, takes advantage of the observation that the AI =1/2 selection rule 
in kaon decays is well reproduced in terms of three basic parameters (the constituent quark 
mass M and the quark and gluon condensates) in terms of which all hadronic matrix elements 
of the AS = 1 lagrangian can be expressed. 

We have used the best fit of the selection rule to constrain the allowed ranges of M, (qq) 
and (GG) and we have fed them in the analysis of e'/e. 

Nonetheless, the error bars on the prediction of e'/e remain large. This is due to two 
conspiring features: 1) the destructive interference between the large hadronic matrix elements 
of Qq and Qg which enhances up to an order of magnitude any related uncertainty in the final 
prediction (this feature is general and does not depend on the specific approach); 2) the fact 
that large quark-condensate values are preferred in fitting the isospin zero K° — > tttt amplitude 
at 0(p 2 ) (which is a model dependent result). 

Whereas little can be done concerning point 1) which makes difficult any theoretical attempt 
to predict e'/e with a precision better than a factor two, an improvement on 2) can be pursued 
within the present approach. 

Two lines of research are in progress. On the one hand, we are extending the analysis 
to 0(p 4 ) in the chiral expansion to gain better precision on the hadronic matrix elements 
and to determine in a self-consistent way the polinomial contributions from the chiral loops; 
preliminary results indicate that the AI = 1/2 rule is reproduced for smaller values of the 
gluon and quark condensates, thus reducing our error bar, in the direction shown by our more 
restrictive estimate. On the other hand, we are studying the AS = 2 sector to determine 
at the same order of accuracy B K and the K L -K S mass difference by including in the latter 
the interference with long-distance contributions that can be self-consistently computed in the 



present approach. 

Whether this program is successfull may better determine how much of the long range 
dynamics of QCD is embedded in the present approach and increase our confidence on the 
predictions of unknown observables. 
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